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Abstract 



d NEMO. 2003^ is a DVD with 91,520 stellar model atmospheres representing a 
5D grid of modified ATLAS9 atmospheres; the purpose of the modifications 
was to include different treatments of convection and higher vertical resolu- 
tion. In addition, for every model fluxes are provided and color indices for 14 
different photometric systems. Because the model grid contained gaps due to 
non-converging models, we developed and applied a set of 4D interpolation 
routines to complete the grid. 

All the data, which will continuously be updated, can be found on the NEMO 
homepage http://cmis.astro.univie.ac.at/nemo/ and are available via 
DVD. 



Introduction 

As described in Heiter et al. (2002, hereafter H02) several new sets of grids of 
model stellar atmospheres were computed with modified versions of the ATLAS9 
code. The individual sets differed from each other, and from previous ones, 
essentially in the physics used for the treatment of the convective energy trans- 
port, in the higher vertical resolution of the atmospheres, and in a finer grid 
in the (Toff, \ogg) plane. The improvements related to resolution enabled the 
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computation of derivatives of color indices, as well as limb darkening coeffi- 
cients and derivatives (Garrido et al. 2001, 2002; Barban et al. 2003), accurate 
enough for pulsation mode identification, and of smooth pre-main sequence 
evolutionary tracks (Montalban et al. 2004). A description of the modifications 
of ATLAS9 with respect to their treatment(s) of convection was also provided 
by H02. While the grids published by Kurucz (1993b (K93b), 1998) and Castelli 
et al. (1997) favor a rather efficient convection in stellar atmospheres due to 
their calibration based on solar central intensities, the new sets of grids were 
computed for convection models which predict a much lower efficiency, at least 
within the photosphere. The motivation behind this has been to allow a choice 
between the two assumptions, as none of the presently implemented local, ho- 
mogeneous models are able to reproduce all observations, whether for the sun 
or for groups of stars spread over wider regions of the HR diagram (see, e.g., 
Smalley et al. 2002). In addition, model atmospheres with lower convective 
efficiency than those in Kurucz (K93b, 1998) and Castelli et al. (1997) distri- 
butions allow recovering a larger number of observations for a given choice of 
convection model and its parameters (H02, D'Antona et al. 2002, D'Antona Sj. 
Montalban 2003, Montalban et al. 2004, Stassun et al. 2004). Schuler et al. 
(2004) used different sets of ATLAS9 grids, including the 72 layer models pre- 
sented here, to derive oxygen abundances in open cluster dwarfs. Their work 
shows that the differences in abundances due to different models are small 
compared to the typical error in the relative oxygen triplet abundances. 

Parameter range 

As discussed in H02, the computations for all values of the grid parameters have 
been finished and we can present the whole 5D grid (Table 1). 

The convection models CM (Canuto & Mazzitelli, 1991) and CGM (Canuto, 
Goldman & Mazzitelli, 1996) and their implementation into ATLAS9 are sum- 
marized in H02. As is shown there as well, the convective flux can be neglected 
for temperatures higher than 8 600 K. Thus, convection has been turned off for 
such models and they are computed only once for models with 288 layers and 
once for models with 72 layers and are included in the CGM subgrids. This ex- 
plains the two different temperature ranges in Table 1 and the total number of 
models of 91,520. We also repeat here that the uppermost layers of our models 
are located at logrRoss = —6.875 and the difference of consecutive layers in 
logTRoss is 0.125 and 0.03125 for models with 72 and 288 layers, respectively.^ 
For all models, the opacity distribution functions calculated by Kurucz (1993a) 
were used. 



The lowermost layer is located at logTRoss = —6.875 + (72 — l)/8 = +2.0 and 
-6.875 -I- (288 - l)/32 = +2.09375 for models with 72 and 288 layers, respectively 
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Parameter 


From 


To 


Increment 

(List of Values) 


Unit 


# 


T'off 


4,000 


10,000 

8, eoo'^ 


200 


K 


31 

24" 


log 5 


2.0 


5.0 


0.2 


{cgs}" 


16 


[M/H] 


-2.0 


+ 1.0 


(-2.0, -1.5, ±1.0, 

±0.5, ±0.3, ±0.2, 

±0.1, 0.0) 


dex 


13 


^micro 





4 


(0, 1, 2, 4) 


kms^^ 


4 


Convec. 
/rcsol. 


CGM(a*=0.09) / 72, MLT(a==0.5) / 72, 
CGM(a*=0.09) / 288, CM / 288 


4 


Total number of model atmospheres 


91,520 



Table 1: Parameter range of the grid 

'Convec' = convection model, 'resol.' = vertical resolution (number of layers) 

" for the non CGM models, see also text 

'' log cm s^^ 



Colors 



Colors and color indices in 14 photometric systems have been calculated with 
the program "colors", which is based on the programs of Kurucz CD-ROM 13 
(http://kurucz.harvard.edu/prograins/CDLORS/) revised and rewritten in 
FortranQO by JN. This program is also available on the NEMO website as well 
as on DVD. The color indices for specific photometric systems have been nor- 
malized using model fluxes and measured indices for standard stars as follows: 
(3 Leo (roff=8 850, log.g=4.16, [M/i7]=0.0, w„ii„o=2) for Walraven VBLUW, 
HD 83373 (9 250, 4.00, +0.0, 2) for the hkl94 index of Stromgren uvbyCa, a 
B star (30 000, 4.00, +0.0, 2) for Vilnius, and Vega (9 550, 3.95, -0.5, 2) for 
Cousins (Cape) VRI, Aa, Geneva, JHKcit, (V)RIJKL, thirteen, and UBV. DDO 
and HST color indices have not been normalized. The standard star colors and 
resulting zero points, as well as the references for the filter transmission curves 
that have been used, are given in Tables 2-4. Colors in the DDO system 
were calculated with data from Lamia (1982, L82). The HST filter curves were 
taken from ftp://ftp.stsci.edu/cdbs/cdbs8/synphot_tables/. The fil- 
ter transmission curves for all systems are shown in Figs. 1-5. The symbols 
represent the values given by K93b or found in the literature, the lines show 
the results of the (curvature weighted) parabolic interpolations from the Ku- 
rucz programs, which are sometimes convoluted with other response functions 
(atmosphere, optics, receiver). 
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H-Beta (Vega) 



Kurucz 1993b, 1998 (. . . /programs/COLORS/) 



H 



13 



2.903 



2.695 



Using additional response functions 'OneP21', 'Air', and 'Reflet'. 
Output "468nm.dat" are not correct Hp values, (see text below) 



Delta a (Vega) 



Kupka et al. (2003), Crawford & Barnes (1970), Maitzen & Vogt (1983) 



-0.009 



0.600 



b-y 



-0.005 



0.500 



a = 92- (.91 + y)/2, Aa = a - ao[(fo - v)] [B - V); [gi - y)] 
Using additional response functions 'OneP21' and 'Reflet'. 



DDO 



L82 p. 55 



35-38 38-41 



41-42 42-45 



0.000 0.000 0.000 0.000 



45-48 



0.000 



No zero point star given. 
Using additional response functions 'OncP21' 



'Air', and 'Rcfict'. 



uvbyCa (Vega, HD83373) 



Crawford & Barnes (1970), Hauck & Mermilliod (1980), 
Stromgrcn (1966) 



u- 



1.411 



0.802 



uq- 



1.411 



0.780 



-y 



0.004 



nil 



0.157 



Cl 



1.089 



0.500 -0.070 -0.058 -0.080 -0.133 



ClO 



1.089 



hk 



0.194 



Model for HD83373 (hk=0.194): 

(Toff, logg, [M/Hl wniicro) = (9 250, 4.0, +0.0, 2) 

Using additional response functions 'OncP21', 'Air', and 'Reflet' 



Table 2: Names, zero point stars, references, standard and correction values of indices 
of the photometric systems. The last line contains remarks concerning the Kurucz 
programs among others. 

It turned out that the standard output of ATLAS9 for the flux files is meshed 
too coarse in A to give accurate results for Hp. So the color output for this 
system is called "468nm.dat" instead of "beta.dat". 
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Walraven (/? Leo) 



L82 p. 76, Lub & Pel (1977) 



V-B 



0.034 



'0.561 



B-U 



0.436 



1.232 



U-W 



0.108 



0.106 



B-L 



0.198 



0.460 



A = 372 nm: /i(K93b) = 0.468, /i(L82) = 0.488 
A = 462 nm: 6i(K93b) = 0.282, 6i(L82) = 0.287 

Model for /3 Leo: (Tcff, logff, [M/H], i^micro) ^ (8 850, 4.16, +0.0, 2) 

Geneva (Vega) 

L82 p. 61 



U-B 



1.505 



0.759 



V-B 



0.959 



0.293 



Bi-B 



0.900 



0.942 



Bo-B 



1.510 



1.443 



Vi-B 



1.662 



1.045 



G-B 



2.168 



1.318 



K93b has some additional response function values. 



GenevaS (Vega) 



Cramer & Maeder (1979), Kiinzli et al. (1997), L82 p. 61, 
North fc Kupka (1997), Rufener fc Nicolet (1988) 



U-B 



1.505 



0.004 



V-B 



0.959 



0.017 



Bi-B 



0.900 



0.017 



Bo-B 



1.510 



0.005 



Vi-B 



1.662 



0.018 



G-B 



2.168 



0.027 



Additional (derived) indices are calculated: 
U-B2, B1-B2, B2-V1, Vi-G, d, A, g, m2, X, Y, Z, pT, pG 

UBV (USA I Vilnius | Buser) (Vega) 

Arp (1961), Azusicnis & Straizys (1969), Buser (1978), 
Buser & Kurucz (1978), L82 p. 48, Matthews & Sandage (1963) 



U-B 



-0.005 



ui-bi 



-0.475 



ui.3-bi.. 



-0.468 



U3-B2 



-0.408 



B-V 



-0.003 



bi-vi 



0.619 



bo-vo 



0.601 



B, 



-vo 



0.601 



Vilnius (7 color) (Vega) 



L82 p. 66 



U-P P-X X-Y Y-Z Z-V V-S T-S 



0.630 



0.304 



0.830 



0.302 



0.290 



0.543 



0.100 



0.414 



0.050 



0.218 



0.140 



0.826 



0.040 



0.267 



Normalization in L82: all indices are 0.000 for 
(Teff, log 5, [M/H], t;„,iero) = (30 000, 4.00, +0.0, 2) 



Table 3: see Table 2 
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HST ('External') 



ftp : //ftp . stsci . edu/cdbs/cdbs8/syiiphot_tables/ 



555 



0.000 



606 



0.000 



814 



0.000 



No zero point star. 
Response functions of cameras 'pcf, 'wf2f', 'wf3f', 'wf4f', and 
'wfpc2_f' used (same color index for one wavelength region). 
Thirteen (Vega) 

L82 p. 68, Johnson & Mitchell (1975) 



33-52 35-52 37-52 40-52 45-52 52-58 52-58' 



0.052 



0.035 



0.043 



0.013 



0.001 



-0.008 



0.001 



-0.202 



0.259 0.090 0.584 0.377 



0.356 0.390 



52-63 52-72 52-80 52-86 52-99 52-110 



-0.011 



0.664 



0.008 



1.103 



0.010 



1.451 



-0.004 



1.683 



0.011 



2.019 



0.000 



2.441 



58'(K93b) = 58(L82), 58(K93b): no reference foimd; 

(52-58)vcga = -0.008 (K93b) vs. 0.001 (L82) 

Cousins (Cape) (Vega) 

Cousins (1976), L82 p. 47, 79, 82 



V-R 



-0.009 



0.593 



V-I 



-0.005 



1.228 



V(K93b)=vo(UBVvii„ius), 

R(K93b)=rc(L82), I(K93b)==ic(L82) (== 0.84 vs. 0.86 at A = 840 nm) 

additional values for RCA31034 ('Cathode75'): 

(A,5'(A)) ^ (880, 0.03), (890, 0.01) 

JHKcit (Vega) 

Kurucz 1993b, 1998 (. . . /programs/COLORS/) 



V-J 



0.030 



19.181 



V-H 



0.030 



16.444 



V-K 



0.030 



14.226 



Using additional response functions 'atmos', 'detec', and 'block' 
(V)RIJKL (Vega) 

L82 p. 71 



V-R V-I V-J V-K V-L 



0.000 



0.717 



0.000 



1.488 



0.000 



2.651 



0.000 



4.850 



K93b uses V of L82 



0.000 



6.699 



Table 4: see Table 2 
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Figure 1: Response functions of four out of the 14 photometric systems (1/4). The 
symbols represent the filter values, the lines show the results of the convolution of 
their curvature weighted parabolic interpolations with other response functions (at- 
mosphere; 'Air' (not in Aa), optics: 'Reflet', receiver: 'OneP21'). 



Missing models - grid gaps 

The grid computations were performed with the perl package SMGT (Stellar 
Model Grid Tool), described in Schmidt (1999, S99). Due to bad initial con- 
ditions or because the parameter region is physically and/or mathematically 
'difficult', 1,110 out of the 91,520 models did not converge with respect to the 
SMGT convergence criteria (S99; H02, their Table 3). Their distribution in the 
grid is shown in Figures 6-10 (p. 11 - 15). For the single parameters we find 
difficult regions for low Toff values (Fig. 6), high (and low) log g values (Fig. 7), 
low \MllI\ values (Fig. 8), and for CGM288 whereas MLT72 is 'easiest' among 
the different convection models (Fig. 10); Wmicro seems not to be as relevant for 
creating problems to achieve convergence (Fig. 9). A more detailed insight — for 
combinations of two parameters — can be obtained from the 3D plots of Figures 
11 - 16 (p. 16 - 21). More figures are available on the NEMO website. 

The total numbers of models fulfilling the primary and secondary SMGT 
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Figure 2: Response functions (filter values and their curvature weighted parabolic 
interpolations) of four out of the 14 photometric systems (2/4). 



convergence criteria are given in Table 5 (p. 23). In this sense the models are 
also called fully and partially converged. Note that before interpolation (next 
section) about 3/4 of the (90,410) models are fully converged and 1/4 partially, 
whereas with the (527) interpolated and consecutively converged files it is the 
other way round: 1/4 fulfil the primary convergence criteria and 3/4 only the 
secondary ones. 



Interpolation 

A Fortran90 program was written to close the gaps in the grid through in- 
terpolation (ipl). Since an extended version of the program should allow the 
interpolation of models (and their fluxes and colors) between the grid nodes 
(Ip2), the following requirements are desired: 

1. Multidimensionality: As there are four numerical parameters — Toff, \ogg, 
[M/H], and Wmicro (the convection model is non-numerical) — a 4D in- 
terpolation is necessary. 
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Figure 3; Response functions (filter values and their curvature weighted parabolic 
interpolations) of four out of the 14 photometric systems (3/4). 
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Figure 4: Response functions (filter values and their curvature weighted parabolic 
interpolations) of two out of the 14 photometric systems (4/4). JHKcit is convoluted 
with the additional response functions 'atmos', 'detec', and 'block'. 



This is realized by means of four ID interpolations, one for each numerical 
parameter, because in the surroundings of a gap there are frequently 
further gaps. Hence, in general it is difficult to get a suitable 4D subgrid 
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Figure 5: Wavelength range of the photometric systems 



neighborhood. 

2. Precision: mmag in colors. 

Tests (comparisons of interpolated values with SMGT calculated ones) have 
shown that linear interpolation results in errors of some ten millimags, so 
at least parabolic interpolation is necessary to shrink down the deviation 
to the desired accuracy of one millimag. 

3. Local procedure: A global procedure would involve 25,000 grid nodes and 
thus make the program complex and slow, and it requires a complete grid 
without gaps. The physics does not demand a global method either. 

This excludes splines which would satisfy the precision and differentiability 
requirements. 

4. Differentiability: First derivatives of color indices with respect to Toff and 
log (7 are used for pulsation mode identification (Garrido et al. 1990); they 
should be computable and continuous at the grid nodes. 

Unlike curvature weighted parabolas as are used in the Fortran77 color 
programs by Kurucz (K93b), distance weiglited parabolas can provide 
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Figure 6; Number of missing models versus Tg-g (added up for the other four parame- 
ters) before and after the interpolation. The two horizontal lines indicate the average 
values per T^s step of 35.8 (=1,110/31) and 18.8 (=583/31). 



continuous derivatives at the grid nodes and fulfill also the other require- 
ments. 



Distance Weighted Parabolas for Interpolation 

Fig. 17 (p. 22) illustrates the interpolation procedure in ID. For an interpolation 
in the interval a;2 < a; < a:;3 we need two left neighboring grid nodes xi and X2 
and two right ones X3 and X4. An interpolating function h(x) is constructed 
for the unknown function f(x) using the four given data points {xi,fi),i = 
1, ... ,4, where 

.n = .nx^) (1) 

We define a (quadratic) 'left' parabola pl{x) by means of 2:1,2,3 and a 'right' 
one pr{x) by means of X2.3a, and for comparison a cubic polynomial pc{x) by 
means of all four; 



Ph{xi) = fi=pc{xi) 

Pl{x2) = f2=Pc{x2) =Pk{x2) 



(2) 
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Figure 7: Number of missing models versus logg (added up for the other four pa- 
rameters) before and after the interpolation. The two horizontal lines indicate the 
average values per logg step of 69.4 (=1,110/16) and 36.4 (=583/16). 



Phixa) 



fi =pc{xi) 



Pk{x3.) 

■ PR (2:4) 



The interpolating function for X2 < x < xt, \s defined as distance weighted 
mean of the two parabolas: 



h{x) 



(3) 



w-LJx) -p-LJx) + wk{x) -pnjx) _ 

Wh + WR 

(x3 - x) ■ pLJx) + {x - X2) ■ Pr{x) 

X3 - X2 

A30 -Phix) + A02 -pnix) 
A32 



which is a polynomial of order 3 in x. Note, that this method for interpolation 
is not restricted to equidistant grids and that wl + wr is constant with respect 
to X. (The definition of A^ is given in Eq. 8 on page 15.) 
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Figure 8: Number of missing models versus [M/H] (added up for the other four 
parameters) before and after the interpolation. The two horizontal lines indicate the 
average values per [M/H] step of 85.4 (=1,110/13) and 44.8 (=583/13). 



It holds by (2) that 

HX2) ^Pl{x2) == ./2 {= Pr{x2) = PC{X2)) 

K^s) ^Pr{x3) = fs {= Phixs) = pcixs)) 
For the first derivative we obtain 

-1 ■ pLJx) + (x3 - x) ■ p'-^jx) + 1 ■ pr{x) + (x - X2) ■ p'yjJx) 



(4) 



h'{x) = 



X3 - X2 

'Wl{x) ■ p'l^ix) + wji{x) ■ p'^{x) [PR-Ph]{x) 



Wh +WR 



X3 - X2 



(5) 



which can be interpreted as the sum of the distance weighted derivatives and 
the interval averaged 'slope' between the right and left parabolas. The second 
term vanishes for x ~ X2 and x = x^ and hence 



h'{x2) = p'l{x2) 

h'{x3.) = PRiXi) 



(6) 
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Figure 9: Number of missing models versus i^micro (added up for the other four 
parameters) before and after the interpolation. The two horizontal lines indicate the 
average values per Umicro step of 278 (=1,110/4) and 146 (=583/4). 



Analogously, if we define left and right parabolas pi,{x; X2,3a) and Pr(x; 
2^3.4,5) (with p\_{x) = Pr{x) because of the same three underlying grid nodes 
X2, X3, and X4), and the interpolating function h{x) as the distance weighted 
mean of the two parabolas for the interval 2:3 < x < X4, we obtain 



h'{x3) 
h'ixi) 



PrM 



PRi^s) 



(7) 



Eq. 7 shows the continuity of the first derivative of the interpolating function 
h'{x) at the grid nodes, and as (5) yields a uniquely defined quadratic poly- 
nomial in-between grid nodes, the continuity of the first derivative as obtained 
from distance weighted parabolic interpolation (2)-(3) holds for the entire of 
xi < X < XN for which interpolation is required. Due to this property and the 
fact of being third order in x, (2)-(3) is sometimes also called a "local spline 
interpolation". 
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Figure 10: Number of missing models versus convection model (added up for the other 
four parameters) before and after the interpolation. The two horizontal lines indicate 
the average values per convection model of 278 (=1,110/4) and 146 (=583/4). 



Using the definitions'^ 

Aoj(x) = 
A„- 



-Aji = Ajfe + Akj = Aoj(xi) 



and Lagrange polynomials [i, j,k — I, ... ,3 or 2, ... , 4),^ 

Aoj Agfe _ X — Xj X — Xfe 



Lijk[x) 



A„- A. 



ij ^ik '^i '^j ^i ^k 



Likj \x) 



1, I ^i 
Ltjkixi) = S^i ^ i' . , I e {ij,k} 



(8) 



(9) 



Pl(x) = fiLi23{x) + .f2L23i{x) + fsLsMx) (10) 

Pr{x) = hL23A{x) + f3L342{x) + UL423{x) 



^We will never use it in the following ambiguous manner: xq — xj = Aqj ^ Auj(xi) ■ 
Xi — Xj , except for the unique case i = 0. 

'^Lijkiix) = (Aoj/Aij) ■ (Aofe/A,fc) • (Aoi/Au), 

Pci^) = /l-'^1234(a;) + f2L234l{x) + /3i3412(x) + f4L4l23ix) 
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Figure 11: Number of missing models versus Toff and logg (added up for the other 
three parameters) before the interpolation. The contour lines are drawn for the half 
and the double of the average value of 2.24 (=1,110/31/16) which itself is indicated 
by the bordering dashed line in the 3D part of the plot. 



we obtain 



h{x) = -: \flWL{x)Li23{x)+f4WK{x)L423{x)- 

^32 '- ^^ V ' ^^ V ' 

Li(x) £4(2:) 

+ /2 [wh{x)L23l{x) + WB.,{x)L234{x)] + 

L2(x) 

+ h [w'l{x)L312{x) + wr.(x)£342(x)]^ I = 

/2A30 



— — 

A32 i A32 



-.fsAc 



U(x) 

A30A01 A40A02 

A21 A42 

A30A01 A40A02 

A31 A43 



(11) 



(12) 
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Figure 12: Number of missing models versus Toff and logg (added up for the other 
three parameters) after the interpolation. The contour lines are drawn for the half 
and the double of the average value of 1.18 (=583/31/16) which itself is indicated 
by the bordering dashed line in the 3D part of the plot. 



-A,nA 



30^02 



A, 



30 



-A31A2 



+ /4 



A, 



02 



A4.SA 



43^42 



where we have taken into account that we also want to interpolate vector 
functions in the model grid, because every fi represents a model atmosphere 
with up to 6-288 — 1728 numerical values. Therefore, the formula is optimized 
for computational performance''' such that each fi occurs only once (and all 
Aij > for xi < a:;2 < a; = xo < xa < a;4). 



Interpolation at grid nodes (Ipl) 

Formula (11) is used for closing the gaps in the model atmosphere grids. Assume 
a missing model at 

^0 = (to, go, zo, Vol Co]) (13) 



■''Fortraii90 array arithmetic 
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Figure 13: Number of missing models versus Tcff and [M/H] (added up for the other 
three parameters) before the interpolation. The contour lines are drawn for the half 
and the double of the average value of 2.75 (=1,110/31/13) which itself is indicated 
by the bordering dashed line in the 3D part of the plot. 



where 



t 


= ^(1) 


— Tcs 


9 


^X(2) 


==log5 


z 


=- X(3) 


= [M/H] 


V 


= X(4) 


^micro 


c 


= ^(5) 


= Convection model resolution 



(14) 



Note, that c = X(5) is a non-numerical parameter of the grid and hence can 
not be used for interpolation. Therefore, we have to consider four 4D subgrids 
(for each value of c) and neglect c hereafter when talking about interpolation. 



For each 



Hs) 



s = l, 



,4, an independent ID interpolation is done in the 



following manner: 

First we look for the next two left and right neighbors in the grid. Their 
parameters and the corresponding (vector) function values® are: 

®In the case of ATLAS9 model atmospheres the quantities representing the atmospheric 
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Figure 14; Number of missing models versus Teff and [M/H] (added up for the other 
three parameters) after the interpolation. The contour lines are drawn for the half 
and the double of the average value of 1.45 (=583/31/13) which itself is indicated 
by the bordering dashed line in the 3D part of the plot. 



for 5 = 1: X(s) 

left : 
right 

f^(l) = 

for s = 2: x^^j 

left : 
right : 

f^(2) = 

for s = 3: X(^\ = 



X{i) =t ih <t2 <to <h < U): 

^i(i) = {ti,9o,za,vo), ^2(1) == {h,ga,zo,vo) 
^3(1) = {t3,9o,za,vo), ^4(1) = (i4,5o,zo,i'o) 
/(^j(i)) == f{U,9o,za,va) 
X{2) = 9 (.91 < .92 < 50 < 53 < 54): 
^1(2) = (io,. 91, 20,^^0), ^2(2) = (io,52,zo,i'o) 

^3(2) = (to, 53, 2^0, I'd), ^4(2) = (io, 54,^0, I'd) 

./(^j(2)) = f{tQ,9t,za,va) 



(s) 



H3) 



= z {zi < Z2 < za < Z3 < Z4): 



left 



^1(3) = {ta,9o,zi,vo), ^2(3) = (io, 50,^2, Wo) 



(15) 



(16) 



(17) 

structure arc the integrated mass, the temperature, the gas pressure, the electron density, 
the Rosseland mean opacity and the radiative acceleration at each depth. 
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Figure 15: Number of missing models versus log p and [M/H] (added up for the other 
three parameters) before the interpolation. The contour lines are drawn for the half 
and the double of the average value of 5.34 (=1,110/16/13) which itself is indicated 
by the bordering dashed line in the 3D part of the plot. 



right : ^3(3) = (io, go, 23,1^0), ^4(3) = (^o,5o, 24,-^0) 
/j(3) = .fiXi{3)) = fito,go,Zi,vo) 

for s — 4: X(s-) = a;(4) = v {vi < V2 < vq < v^ < V4): 

left : xi(4) ^ {to,go,zo,vi), X2(4) ^ (ta,go,zo,V2) (18) 

right : X3(4) = {to,go,Zo,V3), ^4(4) = (to,5o, 20,^4) 
/i(4) = f{xi(4)) ^ f{tQ,go,zo,Vi) 

Subsuming the third lines we can write^ 

ks) = m^s)) (19) 

'^One remark about notation: x = (x(i\,x/2)y^(3)y^(4)[y^(5)]) = (*: 9i ^ji"!, c]), 
^0 = (^(i)o>a;(2)o.a;(3)o,a;(4)o[,a;(5)o]) = {to,9o,zo,vo[,co]); e.g. s=3: ^i(3) = 
{to, go, Zi, vol Co]) = (a:(i)o>a;(2)o>a:(3)i,a;(4)o[,a;(5)o]). I.e. ^(jji is a component of £"^(3). 
Thus, we could define m(i,\{xi(i,\) = rnx,-.^ (Eqs. 21, 26). 
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Figure 16: Number of missing models versus log p and [M/H] (added up for the other 
three parameters) after the interpolation. The contour lines are drawn for the half 
and the double of the average value of 2.80 (=583/16/13) which itself is indicated 
by the bordering dashed line in the 3D part of the plot. 



With Eq. 11 (p. 16, Aij -^ Ay/gj) the ID interpolation for xu^ at xq can 



be done: 



1 4 . _ 
h(s){xo) = ^(s)(io, 50,^0, Wo) = -r '^ h{8)Lt{x(^s)o) (20) 



32(s) 



We introduce weighting factors to allow for neighbors which are not at the 
next grid positions because there are further gaps in between: 



"(-) 



Hs) 



6 

d(s) 
4 

El' 

4=1 



(21) 



where Jn^;, ,. is the grid position number of X(^s)i (number corresponding to 
xi^s)i when counting every grid step from the minimum to the maximum value 
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Figure 17; Example for distance weighted parabolas as interpolating function com- 
pared to a cubic polynomial. h{x) is only used for x — xq £ [x2,X3]. Within this 
interval it is quite close to pc{x) while outside the difference between the two grows 
rapidly. 



of parameter xi^g))-^ 

Thus, d(s-) is the sum of distances in grid step units of the interpolating point 
a;(s)o from the neighbors X(^s)i (or in point of view of 4D: sum of distances of 
xq from Xii^g-f). In the optimal case it is 2 + 1 + 1 + 2 = 6 and the maximum 
weighting factor is 1.0. 

It can happen that only three neighbors are available (2:1 left:right or 1:2 
or 3:0 or 0:3). This is always the case for X(4) = w = Umicro, because there 
are only four different values in total (0, 1, 2, and 4). In this case, a simple 
quadratic (2"*^ order polynomial) interpolation is done and the weighting factor 
is given by: 

1 4 _ 2 

2 ' J( ^ 



w 



(«)[3] 



Hs)m 



Hs)m 



1 
< - 

- 2 



(22) 



'E.g. for s = 1, X 



(s)i 



Hi)i 



: 4 000: m^ 



1 (2 for 4 200, 3 for 4 400, 



..., 31 for 10 000). Note, that we are considering interpolation at grid nodes (missing 
models) and consequently rnx,-.^ is also defined uniquely. 
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Number of models 


Percentage 


Interpolation 


before 


after 


A 


before 


after 


A 


Primary (fully) 


67,674 


67,804 


+ 130 


73.94 


74.09 


+0.14 


Secondary (part.) 


22,736 


23,133 


+397 


24.84 


25.28 


+0.43 


Missing / Interpol. 


1,110 


583 


-527 


1.21 


0.64 


-0.58 


Sum 


91,520 


91,520 


±0 


100.00 


100.00 


±0.00 



Table 5: Numbers and percentages of models fulfilling the primary and secondary 
SMGT convergence criteria before and after the interpolation. Due to rounding the 
percentages may not sum up to 100.00. 



with a factor 1/2 to account for the lower accuracy of this interpolation due to 
one fewer neighbor and 4 as minimal (optimal) distance sum. 

In case of only two neighbors (2:0, 1:1, or 0:2) linear interpolation is applied 
and 



"'^^'^' ^^v -d^r ^isM - 4 



1 1 
< 



(23) 



With only one neighbor (1:0, 0:1) we use identical interpolation (as an initial 
guess) and 



W(s)ll] 



lY 1 



2 / d. 



1 1 
< - 



(s)[i] »a(s)[i] 



If there is no neighbor at all (0:0) we set 

W(s)[0] = 

The general expressions for n = 1, . . . , 4 are: 



V{s)[n] 



ly " (-n3 + 9n2 -14n+12)/6 
2, 



Hs)[n] 



(24) 



(25) 



(26) 



'^(s)H - /, l™a:(,io - rn^(,-,i\ 



Once the ID interpolations in all four dimensions X(^s) (i. 9, z, v) have been 
completed, the final result of grid gap closing interpolation is given by 



1 * 



w 

W 



(27) 
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Figure 18: Histogram of neighbor/interpolation qualities in the gricd with 1,110 miss- 
ing/interpolated models. For relaxation interpolation there are only 13 possible values. 



w 



s=l 



Thus, wi^g-f not only takes into account the distance of neighbors in one dimen- 
sion but also the weighting among the 4 dimensions. Strictly speaking it should 
be written as W(^s)[n(,^ixo)]- 

q{xo) is the interpolation/neighbor quality. The maximum possible value 
of it is 1.0 (100%) provided that in all four dimensions there is an optimal 2:2 
neighbor situation. In our grid q{xo) is never 0.0,'' hence /i(afo) is always well 
defined. 



This procedure (Eq. 13 - 27) is applied to all missing models in each of 
the four convection model subgrids. There are two ways to do so. We call 
the first one direct interpolation: make one run for every missing model as 
described above and take only converged models as neighbors. The second one 



^That would happen in case of no neighbors at all in all four dimensions. The maxi- 
mum value is 87.5% and not 100% because of the maximal three neighbors in v. 
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Neighbors 


Dimension s : a;(g) 


# 


Distance sum 


Weight 


left: right 


l:t 


2:5 


3:z 


4:w 


n 


c^C^)!"! 


W(s)\n] 


3:0, 0:3 


+ 


+ 


+ 


+ 


3 


3+2+1=6 


1/2 . 4/6 = 1/3 


2:1, 1:2 


+ 


+ 


+ 


+ 


3 


2+1+1=4 


1/2 • 4/4 = 1/2 


2:2 


+ 


+ 


+ 


- 


4 


2+1+1+2=6 


6/6= 1 


4 


{1/3, 1/2, l}t + {1/3, 1/2, llg + {1/3, V2, 1}. + {V3, V2}. = 
— /II II l2 l5 90I 9I ol 92 95 ool oil 


q{S,) [%] 


{33i, 37i, 41|, 45|, 50, 54^, 58^, 62i, 66|, 70|, 75, 83^,87^} 



Table 6: Possible neighbor situations and values of distance sum d(s)[„] and weighting 
factor it)(3)[„] for the four parameters ajj^j and resulting 13 possible values of the 
neighbor/interpolation quality q{xo) for relaxation interpolation. 



we call relaxation interpolation: make one run after the other, and take also 
interpolated models as neighbors/" as long as the (values in the) model files 
change. 11 

Fig. 18 shows the distribution of the qualities in the grid before the inter- 
polation. In the case of relaxation interpolation there are no gaps in the grid 
anymore, because also the interpolated files are valid neighbors, hence there 
are only 13 possible values of neighbor quality (Table 6). 

We save the directly interpolated files and the relaxation interpolated ones 
in different directories and apply SMGT to them. Different numbers of models 
converge (fully and partially) in the two directories. As there are more converged 
files in the relaxation directory, which we therefore rate as the 'better' one, we 
replace the interpolated files in our grid with one of the appropriate newly 
converged ones in following order: 

1. Fully converged from relaxation interpolated model 

2. Fully converged from directly interpolated model 

3. Partially converged from relaxation interpolated model 



4. Partially converged from directly interpolated model 



l^Alternatively, we could save the interpolated files from the current run elsewhere and 
replace all the old interpolated models from the last run with the new ones from the 
current run only after the current run has been completed. But as we are not interested 
in the intermediate results and to save time a recently interpolated file replaces the old 
one immediately. 

l^The comparison is done character by character (csh-script) which corresponds to a 
numerical significance of at most 8 digits, depending on the variable to be interpolated. 
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Figure 19: Relaxation progress in loop 1: Number of identically interpolated files in 
consecutive relaxation runs versus number of relaxation run (added up over the four 
convection model subgrids). There are actually stages with decreasing number. 379 
runs are need to converge the procedure for the 1110 missing/interpolated files which 
makes an average slope of about 3 files per run. 



In this w/ay we can replace 427 (115, 12, 267, and 33 from group 1, 2, 3, and 
4, respectively) interpolated files by converged models. 

Fig. 19 shows the relaxation progress of the 1,100 missing models in our 
grid. 379 runs are needed to converge the procedure. 

Now with this partially new grid^^ we can repeat the whole procedure: inter- 
polating directly and relaxationally, replacing interpolated files with converged 
models. After 6 such loops no further progress can be achieved in our grid. A 
total of 527 (=427+62+20+10+6+2+0 in the 7 loops) out of the originally 
1,110 missing models converged, 583 are left interpolated. (They are called 
"missing models after the interpolation" in Figs. 11 - 16.) 

The relaxation progress of the other loops (2 - 7) looks similar to that of 
loop 1 (Fig. 19). 

328 out of the 527 models converged both starting from a relaxationally 



^^ Converged files contain diff'erent values compared to the interpolated ones taken as 
their SMGT starting models. 
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interpolated model and a directly interpolated one, 137 converged only in the 
relaxation case and 62 only in the direct one, so without relaxation interpolation 
only 390 would have converged (720 left missing/interpolated) and without 
direct interpolation 465 (645). 

To sum up, the whole grid gap interpolation consists of seven loops, each 
loop of one direct interpolation and one relaxation interpolation, each relaxation 
interpolation of some hundred runs,^"^ each direct interpolation and each run 
consists of four ID interpolations for every missing model. This gives a total 
of 5,893,792 ID interpolations (of vectors/arrays with 6-288 = 1728 or 6-72 
= 432 components). After each loop SMGT is applied to the interpolated files, 
altogether 4,774 times. 

Interpolation between grid nodes (Ip2) 

The Ip2 is real 4D, it expects a complete 4* (256) neighborhood^"' of grid 
nodes (existing model files, converged or interpolated) 

fijki = lixijki) = I{x(i)i,X(2)],X(ji)k,xi^^)i) = f{U,gj,Zk,vi) (28) 

where i,j,k,l — 1,...,4 around the interpolation point xq — (tQ,gQ,zo,vo) 
which now is not a grid node but still we assume xi^g-^i < a;(s)2 ^ xi^g-^Q < 
X(s)3 < 2;(s)4 for s = 1, . . . , 4, i.e. 2 left and 2 right 'neighbors', in each of the 
4 'directions' (dimensions, parameters) xi^^y 

Using again Lagrange's formalism for polynomial interpolation (cf. Eq. 8, 
p. 15) 

(29) 



^Oj(s) = 


= X(s) - X(s)j 


Aij(s) = 


= X(s)i - X(^s)3 


(-y, .\ - 


Aoj(s) ^i)k(s) 



16 (2"^) left-right combined parabolas can be defined, e.g. (cf. Eq. 9, p. 15): 

PLRRl(^) = fl221 ■ Li23(t) ■ L234{g) ■ L234{z) ■ Li23(v) + (31) 

+ /l222 • Li23{t) ■ £234(3) • i234(^) ' i23l(w) + 

+ /1223 • ^123(0 • -^234(5) • -^234(2;) • L3i2{v) + 

+ /l231 • ^123(0 • ^234(9) ■ L342{z) ■ Li23{v) + 



"379, 365, 340, 331, 280, 216, and 166 in the loops 1 to 7, respectively. 
^■'The base 4 is due to i = 1, . . . , 4 (4 neighbors), the exponent 4 is due to s = 1, 
(4D). 
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/l232 • -^123 (i) • -^234(3 
/1233 • ^123(0 • -^234(5 
/l241 • -^123(0 ' ^234(5 
/l242 • -^123(0 ' -^234(5 
/1243 • Li23{t) ■ -^234 (5 

/344I • -^312(0 ' -^423(5 
/3442 • -^312(0 ■ ^423(5 
/3443 • ^312(0 ■ -^423(5 



)-i342W 


L23l{v) 


)-i342(^) 


L312{v) 


)-L423{z) 


Li23{v) 


)-U23{z) 


L23l{v) 


) ■ L423{z) 


Lzi2{v) 


)-U23{z) 


Li23{v) 


)-U23{z) 


L23l{v) 


) ■ Li2?.{z) 


L^uiv) 



Li 



23 (t) X |i234(5) X ^234(2) X (32) 

iLi23{v) ■ /122I + L23l{v) ■ /1222 + -^312 (f) • /1223J + 

+ ^342(2) X 

iLi23{v) ■ /123I + -^231 (l') • /l232 + -^312 (l') ' /l233J + 

+ £423(2) X 

(-^123(^') • /l241 + L23l{v) ■ /l242 + -^312 (f) ' /l243 , 
+ -^342(3) X L234(z) X 



X (-^123(1') • /l321 + -^231 (^') • /l322 + ^312(1') ' /l323 ) + 

+ £342(2) X 
X (£l23(^') • /l331 + L23l{v) ■ /l332 + £312(1') • /1333J + 

+ £423(2) X 

X (£l23(^') • /l341 + £23l(^') • /l342 + £312 (l') ' /1343J + 
X f £123 (^') • /l441 +£231 (l') • /l 442 +£312(1') • /l 443) 1 + 



+ £312 (t) X |£423(5) X £234(2) X 
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iLi23{v) ■ f 34,41 + L23l{v) ■ ^442 + ^312(1') ' /3443J | 



The complete formula contains 81 (3"*) fijki'- * = 1, 2, 3, j ^ 2,3,4, k = 2, 3,4, 
I = 1,2,3. The order of bracketing [v is innermost) is arbitrary. The distance 
(4D hypervolume) weight for this example is 



WhRRhix) = A30(i) • Ao2(2) • Ao2(3) ' ^30(4) = 

= wi^{t) ■ wnig) ■ wr(z) • wl(w) = 

= (ts - t){9 - 92)(z - Z2)(V3 - v) 

The general expressions are given by 



(33) 



PiYkii^) 
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Anr 



or(4) 



i'¥=i j'^j k'=ik I'^l 



"l'2 



^n[,i,(9) 






Lii'^i!^(v) 



3+i 



3+J- 



= J2 ^-i*i(*) >< { J2 ^J^IJ^(ff) 



i=l+i 



J^i+J 



3+k 



3+1 



Y^ Lkk[k!,{z) X Y^ Lii'^i'^{v) X /ijfe. 



^i5fc[(^) = "^iW • w-Ag) ■ Wf^{z) ■ Wiiv). 



(35) 



where i reaches the value L = in case of a left parabola (in x^ij = t) and 
R = 1 in case of a right one, and analogously for j {x(2) ~ g), k {x[3) = z), 
and I (a:;(4) ~ v), and 



WL(a;(s)) = A3o(s) 
WR(a;(s)) 



A 



02(s) 



= E E E E ^Ji-feK^) = 

i=L,R 3=L,R fc=L,R f=L,R 
4 
= n ^32(s) = (^3 - t2){g3 - g2){z3 " Z2)(l'3 " ^2) 



(36) 
(37) 



s=l 
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For the interpolating function we obtain (cf. Eq. 3, p. 12) 

1111 
W 



M-) - ^my>y£^^M^)-p^M^- (38) 

i=o 3=0 fe=o [=0 
_. 1 1 1 1 

^ H "'J^^) Y^ ^0 (5) Y^ Wf^ (z) Y^ wi{v) ■ Pi-0{t, g, z, v) = 

i=0 3=0 k=0 1=0 

-.1 1 1 1 
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3+1 3+3 3+fc 3+f 
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3+i 
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1 3+fe 1 3+;" 

Xl^fc(^) XI -^fcfci^(^)-5I"'K«) XI ^ll'Ai^) ■ fijkl 



k=0 k=l+k 



1=0 1=1+1 



Hijk(v) 



Hijiz.v) 

This can be transformed into a sum of (tetracubic) separable expressions again 
by the definitions 



L2{x) ^ W-lIx) ■ L2'il{x) + WYi{x) ■ L2S.i{x) 
Lsix) = w-l{x) ■ L3i2{x) + W-[i{x) ■ L342{x) 
Li{x) = W-[i{x) ■ L423{x) 



(39) 



And so 
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= Wl{v) lLi23{v)fijkl + L23l{v)fijk2 + L3i2{v)fijk3\ + 

+ Wu.{v) lL234{v)fijk2 + L342{v)fijk3 + -Z^423(^')/iife4 | = 
= fijklW-L{v)Li23{v) + 

+ fijk2 {w-L{v)L23liv) + Wji{v)L234:iv)} + 
^ ^ ^ 

l2(v) 

+ Iijk3 {w'l{v)L312{v) + Wr(w)£342(w)} + 

^^ V ' 

+ I I jki W-liiv) L i23{v) ^ 

'• -V ' 

Li{v) 

4 

= Y,Li{v)-f,,u (41) 

1=1 



Analogously 



15 



1 3+fe 

Htjiz,v) = J2^ki^) Yl Lkk[k'^{z) ■ Htjkiv) = (42) 

k=0 k=l+k 

4 

fc=l 
1 3+5 

H,{g,z,v) = J2'^j^9) Yl LnU!,i9)-Hrjiz,v) = (43) 

]=o j=i+'j 

4 

1 3+4 

H{t,g,z,v) = ^W'^{t)^ Lu'^e^{t)-H,{g,z,v)^ (44) 



i=0 i=l+i 



^^In Fortran90 the (one line) command for Hiji^(v) — > Hij{z,v) is e.g.: hlj ( : , : , :) 
= Lz(l)*hijk(: , : , : ,1) + Lz(2)*hijk( : , : , : ,2) + Lz(3)*hijk( : , : , : ,3) + 
Lz(4)*hljk(: , : , : ,4), where the first ':' stands for the ^ of H, and the second 
and third for i and j, respectively; no DO loops have to be written for this, fe = f , . . . , 4 
is written explicitly. 
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Hix) = ^Z,(t).i?,(g,z,«) 



i=l 

m - ^H{x) (45) 

Thus, we finally can write (cf. Eq. 11, p. 16) 

4 4 4 4 

'^(^) ^ ^T.llllllf^i^i~^^^^)^MLk{z)Li{v)^ (46) 

i=i j=i k=i 1=1 



" w XI ^* W 51 ^1^9) 51 ^k{z) ^ Li{v)J^jki 



i=l j=l k=l 1 = 1 



Hijkiv) 



H,A^,v) 



H,(g.z,v) 



H(t,g,z,v) 

The separability is of great advantage for the differentiation, e.g. (s — 2): 

^ - ^ttttr.Mtf^L,iz)Uv)^ (47) 
y i=l J=l fe=l 1=1 ^ 



i=l j = l fc=l i=l 



Hijkiv) 



Hij(z,v) 



H\''\g.z,v) 



H(s)(t,g^z,v) 

with 

wiix) = -1 (48) 

ti;i,(a;) = +1 

Z'i(a;) = wi^ix) ■ L^^six) -Li23{x) (50) 
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L'2{x) = w^ix) ■ L23-i^{x) + wn{x) ■ L'234^{x) - L23i{x) + L23iix) 

L'^{x) = WL(a;) -^312(2;) +Wr(x) •i342(x) -i3i2(x) +L342(x) 
L'^{x) = Wr(x) •i423(x) + L423{x) 

Derivatives of order higher than one (which are not continuous at the grid 
nodes!) and mixed derivatives can also be computed easily in this way. 

If 2;(s)o = a; e [a;i,a::2] or a; G [x]y^i,X]\[] (no two neighbors on each 
side) then simple quadratic interpolation is applied: pq{x; xi,X2,X3) or pq{x; 
xn~2,xn-i,xn)- 

For grids with fewer or more (numerical) parameters/dimensions X(^s) the 
interpolation works analogously.^^ 

Discussion and outlook 

Improvement of convergence 

The rate of success of the grid gap interpolation procedure described in the 
previous section depends on the location in the parameter space. By success 
we mean that newly converged models are obtained when using the interpolated 
ones as starting models. From Figs. 6 to 10 one can see that the procedure 
removed most of the isolated gaps for parameters where most of the surrounding 
models are converged (i.e. with high neighbor quality). It was also successful 
for a part of the transition region in the HR diagram where convection switches 
from thin superadiabatic zones as in A stars to deep quasi-adiabatic zones as 
in solar type stars. This region extends along a line from around 5 500 K 
at log.g = 2 to around 7 500 K for log^ = 5 and shifts to slightly higher 
temperatures for increasing metallicity. For the part of this region with high 
logg, the interpolation presumably produces models which are close enough to 
the 'correct' structure so that the subsequent model calculation is safe from 
temporarily converging to a wrong solution. The gap remaining at the low 
log (7 end of this region (for all metallicities considered here) indicates that the 
interpolation is not accurate enough to overcome this 'transition problem' there. 
Interpolation is also ineffective in the corner of our parameter space char- 
acterized by low metallicity, low temperature and high \ogg, where the tem- 
perature correction algorithm applied in ATLAS9 cannot cope with adiabatic 
convection in the atmosphere. The third region where models remain missing 
after the interpolation concerns models with 288 layers and [M/H] ~ +1.0. 
The region is confined to temperatures between 8 400 and 8 800 K and \ogg 



^^The number of neighbors will not change: i = 1, . . . , 4, but the number of dimen- 
sions/'directions'; s = 1, . . . , D. 
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larger than 3.6. For those models the flux derivative errors high in the atmo- 
sphere (logTRoss ^ —5) do not decrease below 10% and continuously increase 
to about 40% in the outermost layers. The flux errors and temperature correc- 
tions are very small (lower than the primary convergence criteria). The same 
phenomenon occurs in the models with 72 layers, but is confined to even higher 
layers (logrposs ^ ^6) and therefore they have been classified as converged 
according to the secondary criteria. We cannot give a definitive explanation for 
this problem at the moment. 

For a part of the models in the grids the gas pressure decreases to values 
below 10"^ dyn cm^^ for the outermost layers. These values are outside of 
the range of the ODF tables by Kurucz (1993a), and ATLAS9 uses the constant 
value corresponding to the limiting pressure instead. For models with [M/H] — 
—2.0 this is the case for a nearly triangular region with Tc« > 6 500 K and 
log<7 < 4.2, where at these limits only the uppermost layer is affected, while 
the extension of this zone increases gradually towards the inner layers when 
moving towards the 'corner' of our parameter space. In the most extreme case 
(Teff, logg) — (10 000, 2.2) it extends to logTRoss ^ —4.5 and the pressure at 
the outermost layer is « 10^"''^ dyn cm~^. The affected region is smaller for 
higher metallicities. This inaccurate description of the opacity at high layers 
nevertheless does not seem to present a problem for the convergence of these 
models, as the temperature corrections, flux and flux derivative errors do not 
show a correlated behavior. ^^ The only exception might be the [M/H] = 
—2 subgrids where the hottest models show increasingly larger temperature 
corrections for low logg (but small flux and flux derivative errors). 

Interpolation routines 

At present the interpolation routines work only with the NEMO grid and ATLAS9 
atmospheres with 72 or 288 layers. Future versions are intended to deal with a 
variety of different grid structures and resolutions in a flexible manner. 
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^^See the figures at the NEMO website http://arns.astro.univic.ac.at/ncmo/dvd/??????/vmicro_?/, 
where '??????' is one of 'cgm072', 'cgm288', 'mlt072', or 'cm_288', and '?' is one of 
0,1,2,4; files 'info_p.ps' and 'info_m.ps'. 
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